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Abstract. We performed a spectroscopic study of the 1809 keV gamma-ray line from A1 decay in the Galaxy using the SPI 
imaging spectrometer with its high-resolution Ge detector camera on the INTEGRAL observatory. We analyzed observations 
of the first two mission years, fitting spectra from all 7130 telescope pointings in narrow energy bins to models of instrumental 
background and the ^®A1 sky. Instrumental background is estimated from independent tracers of cosmic-ray activation. The 
shape of the A1 signal is compared to the instrumental response to extract the width of the celestial line.We detect the ^*A1 line 
at 16 (T significance. The line is broadened only slightly, if at all; we constrain the width to be below 2.8 keV (FWHM, 2cr). 
The average Doppler velocities of ^®A1 at the time of its decay in the interstellar medium (t ~L04 My) therefore are probably 
around 100 km s“*, in agreement with expectations from Galactic rotation and interstellar turbulence. The flux (3.3 (±0.4) 
lO^^'ph cm^^s^'rad^') and spatial distribution of the emission are found consistent with previous observations. The derived 
amount of ^®A1 in the Galaxy is 2.8 (±0.8) Mq. 
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1. Introduction 


The detailed measurement of 1808.65 keV emission from 
Galactic ^^Al is one of the design goals of the INTEGRAL 
mission JWinkler et alJl20()3l). ^^Al gamma-ray s were discov¬ 
ered in 1982 by HEAO-C JlVlahonev et alJloOi . and since then 
have been considered to be direct proof of ongoing nucleosyn¬ 
thesis in the Galaxy. From the 9-year mission of the Compton 
Gamma-Ray Observatory, in particular from the COMPTEL 
imaging Compton telescope, all-sky imaging in t he 1.8 MeV 
gamma-ray line from ^^Al had been obtain ed iDiehLel-al 




Oberlacklll997t Knodlseder et al. 1999t Pliischke et al 


200 ll) . Prom these measurements, we learn that ^®A1 emission 
extends all along the plane of the Galaxy; hence ^®A1 nucle¬ 
osynthesis is a common Galactic phenomenon rather than lo¬ 
cal to the solar system. The irregular structure of the emission, 
alignments of emission maxima with spiral-arm tangent, and 
comparisons with tracers of candidate ^®A1 sources all have 
pointed to the conclusion that massive stars dominate ^®A1 
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The high spectral resolution of Ge detectors of 3 keV 
(PWHM) at the ^®A1 line energy of 1808.65 keV is expected 
to reveal more information about the sources and their location 
through Doppler broadenings and shifts, induced from Galactic 
rotation and from dynamics of the ejected Al as it propagates 
in the interstellar medium around the sources. This should add 
important new astrophysical perspectives, beyond what could 
be learned from the large imaging database of the COMPTEL 
sky survey. 

Significant broadening of the ^®A1 line has been found 
by the GRIS team fr om its Ge detecto r instrument through 
a balloon experiment (iNava et al.lll996l) . The reported intrin¬ 
sic ^®A1 line width of 5.4 keV corresponds to velocities above 
450 km s“* for decaying ^®A1 isotopes, if interpreted dynam¬ 
ically, from Doppler shifts. Various attempts were undertaken 
to understand how su ch high velocities could be maintained o n 
time scales of Myrs jChen et al 111 997t IStiirner & Navalll 999ll . 
It was clear that only rather unusual circumstances such as 
large interstellar cavities, or high fraction of ^®A1 being de¬ 
posited onto grains near the Al sources, could offer an expla- 
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nation. With recent space-based gamma-r ay spectroscopy ex¬ 
periments, in particular frorn the R HESSI jLin et alJ2002l) and 
INTEGRAL iwinkler et alJl2003ll mission s, it became appar- 
may not be valid (l■Smitbll2nn/!l 12004 
The interstellar medium surrounding 
sources of Al may be less extreme; Al is considered an im¬ 
portant probe of these astrophysical sites. 

In this paper, we discuss our analysis of data from the SPI 
spectrometer on INTEGRAL with respect to the intrinsic Al 
line width and to the large-scale properties of Al emission in 
the inner Galaxy, exploiting data from the first two years of the 
mission. 

2. Observations, data and analysis 

2.1. Data selections and standard processing 
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Fig. 1. Exposure of the inner Galaxy for the data as selected 
from orbits 15-259 for our study. At the Galactic center, the 
effective exposure is 4 Ms, at longitudes +50° 0.85 Ms (linear 
scale). 

We analyze data from the SPI instrument, the coded-mask 
imaging spectrometer instrument on INTEGRAL, with its 19- 
element Ge gamma-ray ca mera; the SPI instrument and its per¬ 
forma nce are described in dVedrenne et alJEoOSt iRoaues et al.l 
I 2 OO 3 I) . Eor the results reported here, we use data from the 
INTEGRAL Core Program (which focuses on a survey of the 
inner Galaxy), our Open Program data, and data that were pub¬ 
licly available at the time of this analysis. Our database spans 
INTEGRAL orbits 15-259, and our filtering to obtain data 
free from contamination due to increased-background events 
or other anomalies yields a total observing time of 16.5 Ms 
from this set. Eiltering the data of glitches towards the perigee 
near the radiation belts of the Earth, from solar-flare periods, 
and from otherwise anomalous background conditions was es¬ 
sential in order to obtain a sufficiently homogeneous database 
in terms of background handling. We selected data within or¬ 
bit phases 0.1-0.9 in order to exclude orbital phases near the 
Earth’s radiation belts, we excluded data where the saturated 
Ge count rate exceeds 10^ counts s ' to reject anomalies, and 
we excluded solar-flare periods with a 1-day window triggered 
by the >30 MeV count rate of the GOES proton monitor ex¬ 
ceeding a value of 0.5 Hz. 


Our standard processing includes calibration of the energy 
scale once per orbit by fitting Gaussians to strong instrumental 
lines. This attains a relative pre cision of 0.05 keV a t 1809 keV 
over these 1.5 years of data jEonioii et al.ll2004li . The non- 
symmetric line profiles, which develop from degradation of the 
Ge detectors due to cosmic-ray bombardement, result in shifts 
of the energy scale of a few tenths of keV at 1809 keV. Due 
to calibration-line fitting with symmetric Gaussians, these de¬ 
tector degradations mainly result in slight broadenings of the 
effective instrumental line resolution; the offset in the absolute 
energy scale should be below 0.1 keV. Eor each pointing of the 
instrument, standard processing assembles the energy spectra 
of each of the Ge detectors at 0.5 keV binning, together with 
detector deadtimes and various housekeeping data used for data 
selections and for background modelling. 

In this study, we make use of single-detector hits only, leav¬ 
ing coincident hits of more than one detector aside. Although 
this reduces the overall detection efficiency by about 50 %, our 
background and spectral-response modeling is substantially 
more straightforward for this case. Altogether our database 
covers the Galactic plane (pointings within 30 degrees latitude) 
with 7130 pointings (135470 spectra), equivalent to a total 
deadtime-corrected exposure of 12.86 Ms, with, e.g., 4 Ms at 
the Galactic center and 0.85 Ms at longitudes +50° (see EigQ. 

2.2. Background treatment 

Typically, source signals such as Galactic Al contribute only 
at the percent level, the remainder of the total signal being in¬ 
strumental background. Therefore the treatment of background 
is a key aspect of the achieved sensitivity and suppression of 
systematics, as is the case for all gamma-ray instruments. In 
principle, with SPI and its coded mask we can distinguish “sig¬ 
nal” from “background” by comparing the 19-detector ampli¬ 
tude patterns rj - djj 2 j dj between sky observations and a 
suitable reference (j is detector index, d is counts); instrumen¬ 
tal background should be responsible for a pattern indepen¬ 
dent of instrument pointings. But diffuse sky emission such as 
from Galactic ^®A1 is modulated rather weakly by the coded 
mask between successive pointings with typically 2° offsets, 
unlike for strong point sources. Therefore, the “background ref¬ 
erence” pattern cannot be taken from pointings nearby in ob¬ 
serving time, which is desirable to probe a similar background 
situation. Eor our large database, however, which spans a con¬ 
siderable range in time and sky-pointing directions, we can 
detect the celestial diffuse ^®A1 signal when we use the time- 
averaged data for the background pattern reference (see Eig.|3. 
Independent data are preferable, e.g. when the instrument was 
pointed towards high latitudes, or using detector count rates 
outside the energy range of interest, or using independent trac¬ 
ers of instrumental background. Instrumental-background dif¬ 
ferences between such references and actual data may be ex¬ 
pected, though; some time variation in the background model 
should be allowed for by using additional parameters when fit¬ 
ting measurements. 

Eollowing this general idea, we developed a background 
model for our study; more details will be described in a forth- 


ent that the GRIS result 
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coming paper (Halloin et al., in preparation). We base our 
model on presumed tracers of instrumental background, and 
optimize these to then establish the “background” detector ra¬ 
tio pattern per pointing. 

First, we correlate a larger set of potential background trac¬ 
ers with the count-rate variations in our energy band of inter¬ 
est around the ^®A1 line. We choose the broad band of 1800- 
1820 keV for this correlation, and the mission part where all 
19 detectors of the camera were fully functional. We sum over 
all 19 detectors to improve statistical precision for the detector 
data. This allows us to identify correlations of candidate back¬ 
ground tracers with actual background variations in our cam¬ 
era with high significance. The best-correlated tracer is identi¬ 
fied as the “prime” background tracer. Then we remove the in¬ 
formation of this prime tracer from all other candidate tracers 
through orthogonalization, and again correlate the (now addi¬ 
tional to our prime tracer, or orthogonal-only) components of 
the other candidate tracers with our count rates, to iteratively 
find more “next-best” tracers. No more than three orthogonal- 
ized background-tracer components are needed for adequate 
background modelling of the energy bands around the ^®A1 
line. This first step allows us to identify the necessary set of 
tracers and their decomposition order. We find a hierarchy of 
(1) non-saturated Ge camera count rate, (2) cumulative satu¬ 
rated Ge counts since mission start, and (3) plastic anticoinci¬ 
dence count rate. In a second step, this decomposition scheme 
is applied to our spectroscopic dataset in 0.5 keV energy bins, 
separately for three time ranges to account for the two detector 
failures. This determines for each tracer template the detector 
patterns of instrumental background as fitted to the measured 
data. In this way, our background model is based on actual 
background-tracer data with their high statistical precision and 
high time resolution, yet the detector count ratios within the 
19/18/17-element camera and per energy bin may be different. 
In the final model fitting step (see next Section 2.3), only a few 
(3-10) background parameters are adjusted to allow for grad¬ 
ual drifts of signal-to-background ratios with time, and/or to 
account for normalization differences for the times where one 
or two detector elements of our camera (detectors 2 and 17) 
failed. 

2.3. Spectra determination 

Our measured data consist of spectra per detector and pointing 
(the typical spectral signatures are illustrated in the cumula¬ 
tive spectrum shown in Fig|2j. We derive our spectra of celes¬ 
tial emission by fitting our set of model components for back¬ 
ground and a (set of) model (components) for the spatial distri¬ 
bution of celestial emission to our data, fitting amplitude coef¬ 
ficients independe ntly for fine energy bins in the range 1800- 
1820 keV (see also lStrong et al.l200a for method details). The 
count data in our spectra per pointing and detector are modelled 
as: 

Ns Nb 

^ijk - ^ ^7jk ^ ^ ^ ®b,tBb-ijk 

m,n s t b 

with i,j,k as indices for data space dimensions energy, detector, 
pointing, m,n indices for sky dimensions longitude, latitude. 



Fig. 2. Raw data spectrum in the region around the Al line, il¬ 
lustrating instrumental line and continuum backgrounds (from 
single-detector events (SE) during the first inner-Galaxy sur¬ 
vey (GCDE-1)). Line origins are radioactivities excited by 
cosmic ray bombardements, specifically ^°^Bi (1764.3 keV), 
28A1 (1779.0 keV), and ^^Na (1808.7 keV) with ^^Mn 
(181 0.9 keV) for the composite feature underlying the ^®A1 
line jWeidenspointner et alJl200^ . Eitting these instrumental 
lines we typically find centroids and widths of (1764.4/3.1), 
(1779.0/3.2), and (1809.4/4.2), respectively. The intensity of 
the instrumental feature at 1809.4 keV is typically 12% of the 
instrumental line at 1779 keV. 

and Ns sky model components S, Nb background model com¬ 
ponents B. The skymaps S are convolved into the data space of 
measured spectra for each detector using observation attitude 
information and the instrumental response, which had been de¬ 
termined from Monte Carlo simulati ons and was adjust ed to 
prelaunch calibration measurements jSturner et alJl2003l) . We 
assess the quality of the fit through checks on its global preci¬ 
sion (x^) and on the residuals of the fit to our 135470 measured 
spectra. Typical values per degree-of-freedom range from 
1.5 in the instrumental background line at 1779 keV to 0.8- 
1.2 in the Al line region. Typical statistial errors in fitted flux 
values per bin are 0.12 ph cm^^s^'rad^' keV”'. Systematic un¬ 
certainties can be estimated by increasing errors per data point 
until the fit quality criterion x^ yields a value of 1.0; we de¬ 
rive systematic uncertainty to be 40% of the value for statis¬ 
tical uncertainties. Coefficients © for the background intensity 
(time variation allowed) and for the skymap intensity (constant 
in time) are derived, the latter (0^) comprising the resultant 
spectrum of the signal from the sky (e.g. Eig.|3l. 

The celestial origin of the observed feature near 1809 keV 
has been verified through identical analysis on a dataset from 
high-latitude observations (Eig.|3l, where the signal disappears. 
Moreover, the instrumental feature at 1810 keV (see Eig. |2j 
leaking also into the high-latitude data spectrum in Eig. |3l(top)) 
is clearly offset from the observed signal, and is significantly 
broader, due to its origin as a composite of local radioactivities 
jWeidenspointner'et alj2003ll . 

The instrumental line-shape variation due to degradation of 
detectors was determined during the mission by fitting a spe- 
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Energy [keV] 



Energy [keV] 

Fig. 3. Test for spurious instrumental background leakage into 
spectra for celestial emission: Sky model fitting was applied to 
high-latitude data (above), and to data from the inner Galaxy 
(below). The sky intensity distribution was modeled from the 
COMPTEL ^®A1 skymap, and the pointing pattern of the in¬ 
ner Galaxy was adopted also for the high-latitude data. As a 
simple background model we use globally-averaged data; this 
background model still correlates with the measurements as 
taken along the Galactic plane, which causes non-zero inten¬ 
sities at all energies. Nevertheless, the celestial ^®A1 line is re¬ 
vealed from the increased correlation with the skymap around 
1809 keV for the Galactic-plane data. 

cific spectral response to instrumental lines which captures a 
degradation parameter t (see Fig. 0}. This spectral response 
function consists of a Gaussian shape which characterizes each 
particular detector’s intrinsic resolution, and a one-sided ex¬ 
ponential function extending from the peak of the Gaussian 
towards lower energies, which characterizes the pulse height 
losses due to detector degradation through its width r: 



with Si the amplitude per energy bin £,, Eq the photopeak line 
energy, and cr the instrumental resolution of the Ge detector. 



Fig. 4. Evolution of the detector resolution with time for the 
region around 1800 keV. The width t of the one-sided expo¬ 
nential which we use to describe degradation, averaged over 
all 19 detectors, evolves with time (expressed here in 3-day in¬ 
tervals of the satellite orbits) due to degradation and detector 
annealings. The grey areas represent the linear approximations 
which we adopted for modelling detector degradations, with 
their uncertainties (Kretschmer et al., in preparation). 


By fitting many calibration lines over the SPI energy range, the 
time-variable degradation has been determined (see Fig.|3}. In 
this analysis, we use this spectral-response behaviour to accu¬ 
mulate the expected instrumental line shape, which we use in 
fitting the celestial signal in our Al line spectra and to derive 
constraints on the additional broadening caused by intrinsic ve¬ 
locity variances of decaying ^®A1 (e.g. Fig.|6l see below). We 
plan to take this time-variable spectral response directly into 
account in our instrumental response when convolving the in¬ 
put sky, and then also by our modeling of instrumental features 
of the background (Kretschmer et al., in preparation). 

In fitting non-analytical spectral shapes to our spectra, such 
as our time-integrated instrumental-response line shape con¬ 
volved with intrinsically-broadened ^®A1 emission, gradient- 
driven fit algorithms are inadequate - in particular when we aim 
to determine a parameter like the intrinsic Al width, which is 
convolved with the instrumental line profile before being com¬ 
pared to our flux values per energy bin. The probability distri¬ 
bution in this case is very asymmetric, but we wish to perform 
quantitative statistical analysis of our findings. This is not pos¬ 
sible with tools that inherently assume symmetric and smooth 
probability distributions for parameter fluctuations. We make 
use of the Monte Carlo Markov Chain method, preprocesse d 
by simulated annealing llMetronolis et aljri 9-53tTNe^h993li . 
This method migrates through the plausible parameter space 
as guided by prior knowledge about the parameters, but in a 
random-walk like fashion; for large numbers of samples, this 
assembles the parameter value probability distributions rather 
well also for less idalized cases (see Fig. and next Section), 
allowing us to derive an upper limit for the intrinsic ^®A1 line 
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width. Varying priors within plausible ranges has an insignifi¬ 
cant impact on our obtained numerical value. 

3. Results 



Energy [keV] 

Fig. 5. Spectrum derived from sky model fitting using 
COMPTEL Al MaxEnt image and the “orthogonalized trac¬ 
ers” background model. The fit (solid line) uses the shape of 
the instrumental resolution as it results from cosmic-ray degra¬ 
dation and annealings during the time of our measurement 
(see Eig. |3, convolved with a Gaussian for the instrinsic Al 
width, as fitted by MCMC. The ^®A1 line width constraints are 
highly asymmetric (see Eig.|^. Intensities are quoted in units 
of 10“"^ ph cm^^s^'rad^'keV'*. 

By fitting our set of observations at 0.5 keV wide energy 
binning with the sky intensity distribution of ^®A1 as imaged 
by COMPTEL together with our background model based on 
orthogonalized background tracers, we obtain the spectrum for 
Al emission from the inner Galaxy shown in Eig. m The 26 Al 
line is detected at 16cr significance. The line component above 
the linear component of our spectral fit determines the inten¬ 
sity of observed ^^Al. The total ^^Al gamma-ray flux as de¬ 
termined for the inner Galaxy region (-30° < / < 30° and 
-10° < b < 10°, i.e. the flux in the inner radian) obtained 
from this fit is 3.3 (+0.4) 10“"^ ph cm^^s-ij-^jj-i xhis is con¬ 
sistent with both the value of about 4 10^^_gh_cm_^s_^rad^ 
concl uded from previous measurements dPrantzos and Diehll 
1 19961) and in particular the COMPT EL imaging-analysis value 
of 2.8 (+0.4) 10 ph cm-^s 'rad ' ilPliischkel200 1 \ loherlacki 
Il997l) . When we vary the model sky among plausible models, 
we obtain slightly different flux values, within expectations for 
the respective candidate ^^Al source ma ps. The variance of flux 
values in our sample (see Table fTI and iDiehl et al. 1 (I200.5I) ') is 
1.2 10“^ ph cm^^g-ij-ad^* or 4%; we adopt this as our system¬ 
atic flux uncertainty. 

The detection of the celestial ^^Al line is significant (> 3cr) 
in 6 of the 0.5 keV-wide bins covering the ^^Al line; this allows 
us to derive line shape details for the ^^Al line. 


Table 1. ^^Al line results for different sky disttibution models 


Sky Model 

d/l* 

l2 

dl2 

width’ 

exponential disk'* 

n/a 

2.88 

0.33 

1.09 

young disk^ 

-0.62 

3.22 

0.43 

1.19 

free electrons® 

-1.49 

3.54 

0.36 

1.05 

Ne2001( 180)2 

-2.87 

3.22 

0.38 

1.34 

Ne2001(140)* 

-2.22 

3.06 

0.41 

1.22 

dust lOOpm^ 

-1.14 

4.01 

0.42 

1.17 

COMPTEL 2®A1 '0 

-0.70 

3.04 

0.31 

1.17 


* Fit quality: likelihood ratio difference to first model 
2 in units of ph cm^^j-ij-ad^' 

^ mean of celestial line width (see Fig|^ 
exponential disk , scale radius 4 kpc , scale height 180 pc 
^ young disk from iRohin et all scale height 125 pc 

6 free electron spiral arm model from pulsar dispersion measurements 
iTavlor & Cordesll993h . scale height 180 p c 

’’ free electron model from ICordes & Lazidl2002h . no thick disk, scale 
height 180 pc 

* as 2, but scale height 140 pc 

^ IRAS 100 pm, from skyview iWheelock et '^Il99ll) . after subtrac¬ 
tion of Zodiacal light 
*6 ME map from lPlilschke et 


First, the line centroid is determined at 1808.72 
(+0.19(stat)+0.1(syst)) keV from fitting our model shape 
to the observed feature (see above). This is well within 
the laboratory value for the ^^Al line of 1808.65 (7) keV 
(iFirestone and Ekstr6mll20()3) . The variance of line centroids 
for the different models of Table[2is 0.002 keV. 

Second, the ^^Al line shape rather closely resembles the one 
expected from instrument properties such as intrinsic resolution 
and degradations as experienced between annealings. We fit the 
effective accumulated instrumental line shape, convolved with 
a Gaussian for the celestial ^^Al line broadening, to our spec¬ 
trum. The fitted parameters are the line centroid, the intrinsic 
width of celestial ^^Al, the intensity of the line, and two pa¬ 
rameters for the underlying continuum. The intrinsic line width 
appears to be rather small. But this parameter is not as “well- 
behaved” as the others: The parameter probability distribution 
should be symmettic as for a Gaussian, its width then reflects 
the uncertainty of the fitted parameter value. The probability 
distribution for the 2® Al line width, however, peaks at small val¬ 
ues near 0.2 keV, with a minor reduction towards zero intrinsic 
line width, and gradually decreases towards larger line widths 
(Fig.|^. Formal determination of the mean of this probability 
distribution (to yield a “fitted width value”) and its width (to 
yield the “width uncertainty”) is not appropriate. (Nevertheless, 
we show these results in spectra Figures to at least coarsely de¬ 
scribe the intrinsic line broadening). Therefore, we obtain the 
probability distribution from the Monte Carlo Markov Chain 
parameter fitting; we integrate this distribution over the desired 
probability fractions to obtain upper limits for the intrinsic 2® Al 
line width, see Fig. For a 95% probability equivalent to a 
“2 (t” limit, we thus obtain a limit of 2.8 keV from our SPI 
measurement, while the GRIS value of 5.4 keV has a probabil¬ 
ity of 3 10“^. It turns out that the spectral details (line position 
and width) are practically identical for different adopted sky 
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distributions, while Al fl ux values vary by about 4% among 
different plausible models llDiehl et al. 120051) . 



0 2 4 6 

FWHM of Gaussian [keV] 


Fig. 6. Probability distribution for the Al intrinsic line width 
as fitted in Fig.|3 The formal results (Fig.|3i derived from the 
mean and width of this distribution are inadequate descriptions 
of this highly asymmetric probability distribution. Rather, an 
upper limit of 2.8 keV (2cr) applies, and an intrisically even 
much narrower line is fully consistent with our measurement. 


4. Discussion 


Our new measurement of the inner-Galaxy integrated ^®A1 
emission yields a flux that tr anslates into a Gala ctic steady-state 
mass of 2.8 (±0.8 ) Mq (see IPiehl et ai~ll2005l for details and 
implications). From a com parison to earlier measurements (see 
Prantzos and Diehilll996l) . our flux value is on the low side of 
the average among experiments. This trend for imaging versus 
large fleld-of-view experiments may suggest that either imag¬ 
ing helps to suppress instrumental background leakage, or that 
there is some large-scale extended ^®A1 emission from nearby 
sources. 

As a first step towards spatially resolved spectroscopy, we 
separately determined spectra for the inner region of the Galaxy 
(-40° < / < 40°, -10° < ^7 < 10°) using a smooth exponential- 
disk model with scaling parameters 4 kpc (radius) and 180 pc 
(latitude extent), which we split into three longitude segments 
at -10° and 10°. This ob tains spectra for the different seg¬ 
ments dPiehl et al. 1120051 and referenecs therein). Fitting the 
^®A1 line as above, we obtained differences in line centroids 
compatible with expectations from Galactic rotation: the ^®A1 
line is slightly blue-shifted at negative longitudes and red- 
shifted at positive longitudes. This supports the view that the 
observed ^®A1 emission arises from the inner Galaxy region, 
rather than foreground. This reaffirms the above Galaxy-wide 
steady-state interpretation of the measured ^®A1 flux and de¬ 
termi nation of the total ^®A1 mass in the Galaxy dPiehl et al. I 
20051) . Furthermore, these spatially separated spectra also sup¬ 
port the asymmetry of the ^®A1 emission in the inner Galaxy 
indicated in COMPTEL images dPliischke et al.lEoOl]) . in that 
the 4'^‘ quadrant appears brighter than the P'. 

It is evident that line broadenings of a few keV disagree 
with our data from the inner Galaxy (Figs. |3 and |6j; thus, ve¬ 
locities of decaying ^®A1 isotopes are modest. The GRIS mea¬ 
surement of 5.4 keV line broadening appears inconsistent with 


all other measurements to date and is clearly ruled out by our 
result. 

A line broadening of 0.8 keV corresponds to thermal 
Poppler velocities of 100 km s ', as a reference. In the in¬ 
ner Galaxy, Galactic rotation alone leads to Poppler shifts, 
which are particularly prono unced towards longitude s +30°, 
and range up to 0.25 keV dKretschmer et al.l l2003l) : a line 
broadening of about 1 keV had bee n estimated if integrate d 
over this inner region of the Galaxy dKretschmer et alJl2003l) . 
Are our line-shape constraints consistent with Galactic-rotation 
effects and standard interstellar velocities? 
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Fig. 7. Constraints on the width of the Al line from different 
instruments (Icr error bars). 
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We believe so dPiehl et al. II200-5I) . and consider the mea¬ 
sured line width of ^®A1 from the inner Galaxy to be consis¬ 
tent with Galactic rotation and modest (see below) interstellar- 
medium turbulence around the sources of ^®A1. We thus con¬ 
firm earlier re sults obtained by HEAO-C, RHESS I . and SPI on 
INTE GRAL dMahonev et alJll982l: LSmithll2003t Ipiehlet al.l 
l2004h (see Eig. m This is reassuring: neither the suppression 
of deceleration of Al before decay in large, kpc-si zed cavities 
aroun d its sources nor other exotic explanations dChen et alJ 
Il997h are required to account for a large broadening of the Al 
fine. 

How much line broadening from ISM is plausible? Some 
broadening is expected from the apparent existence of rnajor 
inters tellar cavities in the regions of massive-star clusters dOevI 
Il996l) . and from the ejection kinematics of these presumed 
^®A1 sources, i.e. winds and supernovae . WR wind velocities 
are 1200 km s * dVink & deKoteiiE005l) or higher dPliischkel 


120011 ) . and models for ejection of ^®A1 by c ore-collapse super¬ 
novae predict velocities in the same range dHerant & Woosle^ 
11994 . Therefore it will be interesting to test with SPI on 
INTEGRAL whether we can observe such effects for local¬ 
ized regions of Al emission. Here sources may not yet have 
obtained a long-term equilibrium with their surrounding ISM, 
and ^®A1 may thus be preferentially decaying from its initial 
and fast phase, rather than already being slowed down to nor¬ 
mal ISM velocities in the 10 km s * range. The first SPI re- 
sults of a modestly-broa dened ^®A1 line in the Cygnus region 
dKnodlseder et al.l 2004 . as well as the COMPTEL result of 
the Orion region (where ^®A1 emission appears offset from its 
sou rces in the Or ion OBI association towards the Eridanus cav¬ 
ity jp ie hil l2nn2l) T are hints that ^®A1 streaming occurs at un- 
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usually high velocities and kinematics may be peculiar around 
regions with a rather young (few Myr) population of massive 
stars. 

Our constraints on large-scale integrated broadening of the 
^®A1 line from the inner Galaxy can be interpreted in terms of 
interstellar-medium characteristics, if we adopt the contribu¬ 
tion of Galactic rotation from a model: For 1 ke V broadening 
due to Galactic rotation (iKretschmer et alJ|2003l) . the intrinsic 
width from ISM turbulence would be 1.2 keV to obtain our Icr 
line width limit of 1.6 keV. This corresponds to 300 km s“' 
for a 2cr limit on ISM velocities, leaving typical velocities up 
to ^ 100 km s \ well within the acceptable range. Therefore 
we conclude that, within uncertainties, the average velocities 
of decaying ^®A1 in the Galaxy are probably not in excess of 
typical values for the ISM near massive stars. 
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